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FOREWORD 


This  report  summarizes  the  results  of  an  investigation  into  transonic  high  incidence  unsteady 
aerodynamics.  Transonic  wind  tunnel  tests  were  conducted  for  a  semispan  straked  delta  wing 
model  oscillating  in  pitch  at  high  incidences. 

This  test  was  conducted  under  a  cooperative  program  of  research  between  the  Lockheed  Fort 
Worth  Company  (LFWC),  Fort  Worth,  Texas,  USA  (Formerly  the  Fort  Worth  Division  of  Gener¬ 
al  Dynamics  until  28  February  1993)  and  the  National  Aerospace  Laboratory  (NLR),  Amsterdam, 
The  Netherlands.  The  test  was  conducted  in  ApriFMay  1992.  The  models  and  corresponding 
support  system  were  designed  at  NLR  under  an  earlier  separate  program  with  funding  from  the 
Fort  Worth  Division  of  General  Dynamics  (now  LFWC)  and  NLR.  The  fabrication  of  models, 
test  preparation,  wind  tunnel  test  and  reporting  were  performed  at  NLR  under  a  subcontract  from 
LFWC.  This  work  was  funded  under  Air  Force  contract  F33657-84-C-0247  (CCP  4551)  for  the 
Aeronautical  Systems  Center,  Wright-Patterson  Air  Force  Base,  Ohio.  The  work  was  adminis¬ 
tered  by  Mr.  F.  Zapata  of  the  F-16  SPO  (ASCA^PEF)  and  Mr.  L.  J.  Huttsell  of  the  Flight  Dynam¬ 
ics  Directorate  of  Wright  Laboratory  (WL/FIBG)  Wright-Patterson  Air  Force  Base,  Ohio  under 
work  unit  24010292. 

The  program  manager  was  Dr.  A.  M.  Cunningham,  Jr.  at  LFWC.  The  principal  investigators 
were  Dr.  Cunningham  at  LFWC,  Mr.  R.  G.  den  Boer  for  the  wind  tunnel  test  programs  at  NLR. 
Mr.  den  Boer  was  assisted  by  the  following  NLR  specialists:  C.  D.  G.  Dogger,  E.  G.  M.  Geurts, 
A.  P.  Retel  and  R.  J.  Zwaan.  The  authors  would  like  to  acknowledge  the  following  person  who 
contributed  substantially  in  the  project: 

Mr.  E.  W.  M.  Slijkerman  for  the  design  of  the  wind  tunnel  model  and  support;  Messrs.  T. 
Horsman  and  F.  Hofman  for  the  fabrication  and  instrumentaion  of  the  model;  and  Messrs.  O.  van 
Teunenbroek  and  A.  van  der  Kamp  for  their  contribution  in  the  software  development. 

The  test  program  for  the  straked  delta  wing  is  documented  in  three  separate  reports.  This  re¬ 
port  (Part  1)  contains  a  description  of  the  model,  test  setup,  data  acquisition,  and  data  processing. 
Part  2  contains  selected  test  points  for  harmonic  oscillation  and  Part  3  contains  selected  data  for 
simulated  maneuvers.  Parts  2  and  3  are  published  as  WL-TR-94-3095  and  WL-TR-94-3096,  re¬ 
spectively.  An  overview  of  this  test  program  is  also  published  as  WL-TR-94-3017. 
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NONDIMENSIONALIZATION,  SYMBOLS  AND  DEFINITIONS 


For  a  description  of  the  axis  system  and  sign  conventions,  the  reader  is 
referred  to  Appendix  A.  Pressures  and  loads  have  been  made  dimensionless 
according  to  the  following  conventions : 


Nondimensionalization 

Mean  (steady;  the  postscript 

Pressure  :  Cp  0  = 

Balance  loads  :  CN  0  = 

Cm  0  = 

Cl  0 

Sectional  loads :  CN_u  0  = 

CN_1  0  = 

CN_t  0  = 

Cm_u  0  = 

Cm_l  0  = 

Cm_t  0  = 

CN_u  0  = 

Cl_u  0  = 


0  indicates  the  0-th  harmonic  component) 

(p  0  -  ps)  /  q 

Normal  Force  /q*Sref 
Pitching  Moment/g*Sref *cref 
Rolling  Moment/q*Sref *bref 


-  J  cp-  0 

0 

1 

J  cp-0 

0 

1 

I  (Cp-  0  -  Cp*  0) 

0 

1 

-  J  Cp*  0 

0 

I 

J  cp-o 

0 

1 

J  (Cp-  0  -  Cp*  0) 

0 

1 

-  J  Cp*  0 

0 

1 

-  J  cp-  0 


d(x/c)  chordwise 

d(x/c) 

d (x/c) 


(x/c-0.25)  d(x/c) 

(x/c-0.25)  d(x/c) 

(x/c-0.25)  d(x/c) 

d(y/b)  spanwise 

(y/b  )  d(y/b) 


ix 


Unsteady  (all  unsteady  signals  have  been  decomposed  into  haimionic  components; 
the  harmonic  component  is  indicated  by  the  postscript  h; 
each  harmonic  component  has  been  decomposed  into  a  real  (in  phase) 
and  an  imaginary  (out  of  phase)  part;  e.g.) 

Cp  h  =  Re(Cp  h)  +  i*Im(Cp  h) 

Pressure  :  Cp  h  =  (p  h  )  /  q*dalpha 


Balance  loads  :  CN  h  =  Noirmal  Force  /g*Sref *dalpha 

Cm  h  =  Pitching  Moment/q*Sref *cref *dalpha 

Cl  h  =  Rolling  Moment/q*Sref *bref *dalpha 
1 

Sectional  loads:  CN_u  h  =  -  J  Cp*  h  d(x/c)  chordwise 

0 

1 

CN_1  h  =  J  Cp'  h  d(x/c) 

0 

1 

CN_t  h  =  J  (Cp- h  -  Cp*  h)  d(x/c) 

0 

1 

Cm_u  h  =  -  J  Cp*  h  {x/c-0.25)  d(x/c) 

0 

1 

Cm_l  h  =  J  Cp- h  (x/c-0.25)  d(x/c) 

0 

1 

Cm_t  h  =  J  (Cp- h  -  Cp*  h)  (x/c-0.25)  d(x/c) 

0 

1 

CN_u  h  =  -  J  Cp*  h  d(y/b)  spanwise 

0 


Cl_u  h  =  -  J  Cp*  h  (y/b  )  d(y/b) 

0 


X 


Notes : 


Symbols 

acc 

ADC 

alpha, 


1  All  harmonic  (h>0)  components  have  been  nondimens ionali zed  by  the 
first  harmonic  of  dalpha  (in  radians) . 

2  For  layout  reasons,  the  0  indicating  the  0-th  harmonic  component- 
(mean  value)  is  sometimes  omitted  in  tables  and  plots. 

3  Reference  lengths  are  summarized  below: 

Sref  =  0.144406  (m^) 

cref  =  0.820700  (m) 
bref  =  0.417900  (m) . 

4  Pitching  moments: 

-  Wing  :  about  the  rotation  axis . 

-  Sections:  about  25  %  local  chord. 

6  Coefficients  of  spanwise  sections: 

Integration  from  y=0  to  tip;  rolling  moment  about  y=0. 

7  The  section  number  of  the  section  coefficients  is  either  indicated 
at  the  left  hand  side  of  the  presented  values  (e.g.  see  Table  6) 
or  by  an  additional  subscript  (e.g.  CN_u_l) . 

and  definitions 


(m/s=) 


a  (deg) 


(m) 


Acceleration  (acc_ll  is  the  acceleration  measured  by 
accelerometer  11;  see  Table  1) 

Analog  Digital  Converter 

Incidence  relative  to  the  x-axis  as  determined  from  the 
LVDT  signal 

harmonic  oscillations:  zeroth  harmonic  component  of  the 
signal 

maneuvers:  half  the  sum  of  the  maximim  of  the  (1-cos) 
input  and  the  minimum  of  the  cosine  input  of  one 
sequence 

(Note  that  the  incidence  relative  to  the  root  chord  is 
alpha  +  0.0803;  see  Appendix  A) . 

(local)  span;  measured  from  root  chord  of  SiS  model 
{y=0) 


bref 


(m) 


c 

Cd 

Cl 

Cm 

CN 

COP 

Cp 

cr,cref  (m) 


DAC 

dalpha,da  (deg, rad) 


DAT 

dpn,  DPN 
DPT 


f , freq 

(Hz) 

harm,  h 

i 

Im 

k 

(-) 

LCO 

LE 

LVDT 

M , Mach 

(-) 

MCCU 

MCFE 

P 

(Pa) 

p.a. 

P_d 

(deg 

PHARAO 

PLTDA 

reference  span  used  in  nondimensionalizing  the  rolling 
moment:  0.4179  (distance  between  y=0  and  tip  section, 
excluding  the  wing  tip  fairing  (see  Sections  2.2.2  and 
2.2.3) 
chord 

drag  coefficient 

rolling  moment  coefficient 

pitching  moment  coefficient 

normal  force  coefficient 

Computer  of  Pharao  measurement  system 

pressure  coefficient 

length  of  reference  chord 

SiS:  root  chord  (at  y=0)  (  0.820700) 

Digital  Analog  Converter 

model  amplitude  as  determined  from  the  LVDT  transducer 
signal 

harmonic  oscillations:  first  harmonic  component  of  the 
signal 

maneuvers:  half  the  difference  between  the  maximum  of 
the  (1-cos)  input  and  the  minimiim  of  the  cosine 
input  of  one  sequence  (see  also  Figure  12) 

Data  Acquisition  computer  of  Tunnel 
data  point  number 
Data  Processing  computer  of  Tunnel 
frequency,  frequency  of  model  oscillation 
harmonic  component; 

harm  =  1  refers  to  the  model  frequency 

Imaginary  part  e.g.  CN  h=  Re  (CN  h)  +  i*Im(CN  h) 
reduced  frequency,  Jt*f*cref/V 
Limit  Cycle  Oscillation 
Leading  Edge 

Linear  Variable  Differential  Transducer,  refers  to  the 
displacement  transducer  mounted  between  a  fixture  on  the 
turntable  and  a  crank  on  the  main  axis 
freestream  Mach  niomber 

Multi  Channel  Conditioning  Unit  (see  Section  2.3) 

Multi  Channel  Front  End  (see  Section  2.3) 

pressure 

pitching  axis 

pitch  deflection  of  balance  (>0  nose  up) 

Processor  for  Harmonic  And  Random  Oscillations 
Phase  Locked  Time  Domain  Averaging 
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ps 

(Pa) 

q,Q 

(Pa) 

R 

r_d 

(deg) 

Re 

(-) 

Re 

sis 

SITES 

Sref 

(m^) 

T 

(s) 

TE 

V 

WRP 

X 

(mm) 

x/c 

(-,%) 

y 

(mm) 

y/b 

(-,%) 

z 

(mm) 

Greek 

a  (deg) 


V  (mVs) 

Subscripts 


freestream  static  pressure 
freestream  dynamic  pressure 
radius  (see  Figure  6) 

roll  deflection  of  balance  (>0  portside  down) 

Reynolds  number  Re=V*cref/u 

Real  part,  e.g.  CN  h=  Re  (CN  h)  +  i*Im{CN  h) 

Simple  Strake 

Shock  Induced  Trailing  Edge  Separation 
wing  reference  area:  total  wing  area, including  strake, 
excluding  the  area  for  y>0,  see  Figure  6  and  Appendix  A 
(0.144406) 

Duration  of  a  full  (1-cos)  input  (see  Figure  12), 

T=l/3 . 8 
Trailing  Edge 
freestream  velocity 
Wing  Reference  Plane 

ordinate  (see  Appendix  A;  see  also  x/c) 
relative  chordwise  position 

spanwise  ordinate  (see  Appendix  A;  see  also  y/b) 
relative  spanwise  position 
ordinate  (see  Appendix  A) 


incidence  relative  to  the  x-axis 

harmonic  oscillations :  zeroth  harmonic  component  of  the 
signal 

maneuvers:  half  the  sum  of  the  maximum  of  the  (1-cos) 
input  and  the  minimum  of  the  cosine  input  of  one 
sequence 

(Note  that  the  incidence  relative  to  the  root  chord  of 
the  sis  configuration  is  alpha  +  0.0803;  see  Appendix  A) 
(freestream)  kinematic  viscosity 


ELL 

ETD 

m 

c 

s 


Elliptic  (see  Section  4.1  and  Table  3) 

Equal  Time  Delay  (see  Section  4.1  and  Table  3) 
model  (f^  =  model  frequency) 
cut-off  (fc  =  cut-off  frequency) 
sample  (fj  =  sample  frequency) 
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Superscripts 


upper 

lower 


Postcripts 


following  order  postscripts  for  section  coefficients: 
Cl_2_3  h  1 :  N  Normal  force 

m  pitching  moment 
1  rolling  moment 
2 :  u  upper 
1  lower 
t  total 

3 :  section  number 
h  harmonic 


h  harmonic;  when  no  harmonic  is  indicated  the  mean  value  is  presented 

i  instationary 

_in  inertia  part 

_j  section  number 

_1  lower 

_m  mean  (0-th  harmonic) 

_u  upper 

_t  total 

Zero  (O-th  harmonic)  mean;  when  no  harmonic  is  indicated  the  mean  value 

is  presented 

0  (0-th  harmonic)  mean;  when  no  harmonic  is  indicated  the  mean  value 

is  presented 


1  INTRODUCTION 


1 ■ 1  Motive 

The  unsteady  transonic  flow  during  maneuvers  of  fighters  is  not  very  well 
understood.  For  instance,  large  time  delays  and  severe  dynamic  overshoots  in 
normal  force  may  occur,  which  cannot  be  predicted  accurately  by  numerical 
methods.  As  a  consequence  to  be  conservative,  structures  must  be  overdesigned 
or  flight  envelopes  must  be  unnecessarily  restricted.  Therefore,  a  better 
understanding  of  the  unsteady  transonic  vortex  flows  which  occur  during 
maneuvers,  is  of  interest  for  development  and  operation  of  fighters. 

In  response  to  the  above  needs  an  unsteady  transonic  wind  tunnel  test 
simulating  maneuver  aerodynamics  was  conducted  under  a  cooperative  program  of 
the  Lockheed  Forth  Worth  Company  (LFWC,  formerly  the  Forth  Worth  Division  of 
General  Dynamics)  and  the  National  Aerospace  Laboratory  (NLR) . 

This  program  was  funded  by  the  U.S.  Air  Force,  LFWC  and  NLR. 

1 . 2  Objectives 

The  objectives  of  the  unsteady  wind  tunnel  test  devoted  to  maneuver 
aerodynamics  are: 

1  Development  of  a  better  understanding  of  the  physics  of  the  unsteady 
Uiransonic  vortex  flow  about  a  simple  straked  delta  wing. 

2  Generation  of  a  steady  and  unsteady  airloads  data  base  for  use  in  the 
validation  of  fluid  dynamics  computer  codes. 

To  meet  these  objectives,  a  wind  tunnel  experiment  on  a  highly  instrumented 
semispan  simple  straked  wing  model  was  defined.  Both  pitch  oscillation  and 
maneuver  simulation  were  included.  In  a  way  the  test  was  considered  as  a 
continuation  of  the  test  with  a  full-span  model  at  low  speeds  performed  by  NLR 
in  1986  (Ref.  1).  Accordingly,  the  planform  of  the  SiS  model  was  chosen 
indentical  to  the  planform  of  the  full-span  model. 

1.3  Outline  of  this  report 


This  report  deals  with  the  unsteady  transonic  wind  tunnel  test  on  the  semispan 
model  of  the  simple  straked  delta  wing.  Often  it  is  referred  to  as  SiS 
configuration.  The  outline  of  the  report,  which  is  also  indicated  in  the 
foreword,  is  repeated  below. 
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This  report  consists  of  three  parts.  Because  presentation  of  all  the  measured 
data  on  paper  would  require  more  than  30,000  pages,  only  a  limited  number  of 
data  points  is  presented. 

Part  1  presents  the  objectives  of  the  program,  a  description  of  the  model  and 
test  setup,  the  data  acquisition  and  data  processing  techniques  and  also  the 
test  program  for  the  experiment. 

Part  2  presents  five  harmonics  of  five  selected  data  points  of  the  test  on  the 
straked  delta  wing,  as  measured  during  harmonic  oscillation  of  the  model,  in 
plotted  and  tabulated  form.  These  data  are  also  made  available  on  a  floppy 
disk.  As  a  reference,  plots  of  overall  quantities,  deduced  from  a  large 
number  of  data  points,  are  presented. 

Part  3  presents  selected  data  obtained  during  simulation  of  maneuvers  with  the 
semispan  straked  delta  wing.  These  data  are  also  made  available  on  floppy 
disk. 

Readers  who  are  interested  to  obtain  part  or  all  of  the  measured  data  and/or 
the  geometric  data  of  the  model  (e.g.,  in  the  CATIA  or  VDAFS  format)  should 
contact  Lockheed  or  NLR. 
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2  TEST  SETUP 

2 . 1  Wind  tunnel 

The  test  was  performed  in  the  1.60*2.00  High  Speed  Tunnel  (HST)  at  NLR  in 
Amsterdam.  This  wind  tunnel  has  a  closed  circuit  with  a  test  section  length  of 
about  2.5m.  Top  and  bottom  of  the  test  section  are  slotted  walls  with  an  open 

ratio  of  12  percent.  The  velocity  range  of  this  tunnel  is  0<M<1.35  and  by 

changing  the  stagnation  pressure  from  12  kPa  to  390  kPa  a  wide  range  of 
Reynolds  numbers  can  be  covered.  The  diplacement  thickness  of  the  wind  tunnel 
sidewall  boundary  layer  is  7  mm.  For  further  details  the  reader  is  referred  to 
Reference  2 .  Data  in  this  section  refer  to  the  state  of  the  wind  tunnel  during 

the  test.  After  the  test,  in  the  summer  of  1992,  the  tunnel  was  modified.  One 

of  the  modifications  concerned  the  enlargement  of  the  test  section  to 
1.80*2.00  m^,  although  adjustment  of  the  upper  and  lower  walls  still  enable 
operation  in  the  "old"  configuration. 

2.2  Model 

2.2.1  General  arrangement 

For  the  test  on  the  SiS  model,  the  same  highly  instrumented  outboard  wing 
panel  and  support  system  were  used  as  for  a  test  devoted  to  Limit  Cycle 
Oscillations  (LCO) .  Figure  1  shows  the  relation  between  the  SiS  configuration 
and  the  LCO  configuration. 

Figure  2  shows  the  basic  wing  panel  and  corresponding  support  system.  Main 
geometry  information  of  the  wing  panel  is  collected  in  Appendix  B.  Support  was 
provided  through  a  semispan  balance  beam,  mounted  on  a  shaft,  which  in  turn 
was  supported  by  bearings  mounted  on  the  sidewall  turntable. 

Oscillatory  motion  of  the  balance  beam  and  consequently  of  the  outboard  wing 
panel  was  provided  by  an  electro-hydraulic  shaker  system  mounted  on  the 
sidewall  turntable  of  the  tunnel.  Adjustment  of  the  mean  angle  of  attack  was 
done  by  rotating  the  turntable. 

The  sis  configuration  included  a  strake  which  was  attached  to  the  outboard 
wing  (Figure  3) .  In  the  LCO  configuration  a  fighter  type  fuselage  was  fixed 
to  the  turntable. 

All  model  and  support  parts  were  designed  and  fabricated  by  NLR.  Also 
instrumentation  and  calibration  were  accomplished  at  NLR.  For  details  of  the 
design  of  the  two  wind  tunnel  models  the  reader  is  referred  to  Reference  3 . 


turntable 
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Fig.  1  Two  wind  tunnel  models 


bearings 


pitching  axis 


Fig.  3  SiS  configuration  (dimensions  in  mm) 


6 


2.2.2  Outboard  wing  panel,  support  and  excitation 


The  panel  was  linearly  lofted  between  root  and  tip  and  had  a  constant  NACA 
64A204  (mod.)  section.  The  wing  twist  was  -3.0  deg.  Support  was  provided 
through  a  semispan  balance  beam,  moxinted  on  a  shaft,  which  in  turn  was 
supported  by  bearings  mounted  at  the  sidewall  turntable  (Figure  2) .  To  reduce 
inertia  loads,  the  wing  panel  was  made  of  high  strength  aluminum  alloy.  For 
strength  reasons  the  attachment  of  the  wing  panel  to  the  balance  beam  was  not 
at  the  root  chord,  but  at  about  40  per  cent  of  its  span.  The  lower  side  of  the 
wing  panel  was  thickened  over  the  first  part  of  the  span  to  accommodate,  but 
not  make  contact  with,  the  balance  beam.  Main  geometry  information  of  the  wing 
panel  is  collected  in  Appendix  B. 

During  the  design  of  the  model,  the  thickness  of  the  balance  beam  and 
consequently  the  thickening  of  the  wing  were  varied  to  investigate  the 
influence  on  both  the  natural  frequencies  of  the  model -balance  combination  and 
the  aerodynamics.  The  goal  during  this  design  process  was  to  find  an  optimum 
between  conflicting  design  objectives: 

-  A  first  natural  frequency  of  the  wing  panel  and  balance  combination  above 
the  highest  test  frequency. 

-  Minimization  of  the  balance  beam  dimensions  in  order  to  have  the  beam 
stored  within  the  available  interior  space  of  the  wing  panel. 

-  Minimization  of  the  static  deformations  of  the  balance  beam  in  order  to 
avoid  changes  in  the  aerodynamics. 

-  Meeting  the  specified  strength  requirements. 

By  modifying  only  the  lower  part  of  the  wing  in  the  region  between  y  =  -62.3 
mm  and  y  =  -209.9  mm,  a  high  first  natural  frequency  was  obtained,  while  the 
aerodynamics  in  the  outboard  region  were  almost  unchanged.  The  modification 
started  downstream  of  the  leading  edge  flap  and  extended  to  the  trailing  edge. 
The  additional  thickness  varied  from  zero  at  y  =  -209.9  mm  to  2.5  per  cent 
local  chord  of  the  outboard  wing  panel  at  y  =  -62.3  mm. 

The  chordwise  position  of  the  balance  beam  relative  to  the  model  was  chosen  in 
such  a  way  that  the  location  of  the  pitching  axis  was  the  same  as  in  the 
corresponding  low-speed  model  (Ref.  1):  73.3  %  of  the  root  chord  of  the 
semi span  model. 

Oscillatory  motion  of  the  balance  beam  and  consequently  of  the  outboard  wing 
panel  was  provided  by  an  electro-hydraulic  shaker  system  which  consisted  of  a 
hydraulic  power  supply,  linear  actuator  and  servovalve  and  a  feedback  control 
unit  (Ref.  4) .  The  hydraulic  actuator  could  deliver  a  maxim\im  static  force  of 
13000  N  and  a  dynamic  force  of  8000  N  at  a  total  piston  stroke  of  35  mm  at  low 
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frequencies  to  1.1  mm  at  90  Hz.  The  hydraulic  actuator  was  suspended  in  a  box 
which  was  bolted  rigidly  to  the  turntable.  The  piston  was  connected  to  a  crank 
which  converted  the  linear  motion  into  a  rotational  motion.  Because  of 
different  amplitude- frequency  requirements  for  the  LCO  test  and  the  SiS  test, 
the  distance  between  the  working  line  of  the  piston  and  rotation  axis  was 
chosen  80  mm  and  60  mm,  respectively.  The  mean  position  and  dynamic  motion  of 
the  balance  beam  with  respect  to  the  turntable  were  measured  by  an  LVDT 
mounted  between  a  crank  on  the  shaft  and  a  fixture  on  the  turntable. 

The  wing  panel  was  of  a  "clamshell"  design  so  that  all  instrumentation  inside 
the  wing  was  accessible.  Figure  4  shows  the  outboard  wing  panel  with  and 
without  the  balance  beam.  It  had  instrumented  leading  and  trailing  edge 
flaps.  At  the  tip,  a  fairing  with  a  tip  radius  of  0.5  *  the  local  wing 
thickness  was  mounted. 

2.2.3  sis  configuration 

The  sis  configuration  is  shown  in  Figures  1,3,5, 6  and  7.  Since  the  test  on 
this  configuration  was  the  transonic  counterpart  to  the  low-speed  test  of 
Reference  1,  the  strake  had  to  move  with  the  wing  and  was,  therefore,  attached 
to  the  basic  wing  panel.  In  this  way  a  semispan  model  was  obtained  with  the 
same  planform  as  (half  of  the)  full-span  model  used  in  the  test  reported  in 
Reference  1.  Loads  from  both  the  wing  and  strake  were  carried  through  the 
semispan  balance  beam.  The  rotation  axis  was  located  at  the  same  position  as 
in  the  low-speed  test:  73.3  per  cent  root  chord.  The  spanwise  cross-section  of 
the  strake  has  (half)  a  diamond  shape  with  half  a  top  angle  of  11.4  deg.  At 
the  tunnel  wall  the  cross-section  was  rounded  with  a  radius  which  varied 
linearly  from  0  mm  at  the  apex  to  100  mm  at  a  position  417.2  mm  downstream  of 
the  apex.  The  sharp  leading  edge  was  in  the  wing  reference  plane.  The  y=0 
plane,  which  corresponded  to  the  plane  of  symmetry  of  the  (full-span) 
low-speed  model,  was  located  on  a  distance  from  the  sidewall  which 
corresponded  to  the  local  displacement  thickness  of  the  tunnel  sidewall 
boundary  layer  (7.0  mm).  To  impose  the  start  of  the  vortex  on  the  apex,  a 
little  flat  plate  of  1.5  mm  thickness  was  attached  to  the  model  (Figures  5  and 
6) .  Downstream  of  this  plate  the  model  was  extended  from  the  symmetry  plane  in 
starboard  direction  to  y=1.5  mm,  leaving  a  gap  of  5.5  mm  between  sidewall  and 
model.  Downstream  of  x=-96.4  mm,  the  model  was  extended  another  5  mm  to  y=6.5 
mm,  providing  passage  of  the  balance  beam  (Figure  7).  Between  x=-96.4  mm  and 
-520.6  mm  a  filler  plate  of  5.0  mm  thickness  could  be  mounted,  reducing  the 
gap  between  model  and  tunnel  wall  to  0.5  mm.  Close  to  the  kink  in  the  leading 
edge,  the  strake  leading  edge  was  thickened  to  provide  a  smooth  connection 
with  the  outboard  wing  panel.  Also  the  region  where  the  outboard  wing  panel 
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Fig.  4a  The  outboard  wing  panel  without  balance  beam 


Fig.  4b  The  outboard  wing  panel  with  balance  beam 
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Fig.  5b  The  SiS  configuration  in  the  wind  tunnel;  detail  of  apex  with  "little  flat  plate" 


Fig.  5a  The  SiS  configuration  in  the  wind  tunnel;  overall  view 
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Fig.  7  Positioning  of  little  flat  plat  and  filter  plate  (SIS  configuration;  dimensions  in  mm) 
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was  attached  to  the  inboard  part  of  the  model  was  smoothed.  The  leading  and 
trailing  edge  flaps  were  fixed  at  zero  deg.  The  strake  was  provided  with 
three  accelerometers  and  an  additional  spanwise  pressure  section  at  40.42  per 
cent  of  the  root  chord;  for  details  of  the  instrumentation  the  reader  is 
referred  to  Section  2.2.4.  The  model  was  provided  with  transition  strips  of  2 
mm  width  on  the  upper  and  lower  side  of  the  outboard  wing  (y  <-108.65  mm), 
starting  14.5  mm  downstream  of  the  leading  edge  (measured  perpendicular  to  the 
leading  edge)  .  The  grit  size  was  88  vuci  (Carborundiam  150)  . 

Directly  after  the  pressure  and  loads  measurements  on  the  SiS  model,  a 
preparatory  visualization  test  was  done  during  half  a  day  using  the  vapor 
screen  technique  and  a  high  speed  video  camera.  For  this  test  the  large 
markers,  shown  in  Figure  5,  were  painted  on  the  model.  The  results  of  this 
visualization  test  are  not  included  in  this  report;  additional  flow 
visualization  test  are  planned  for  1995. 

2.2.4  Instrumentation 

The  main  wing  balance  measured  normal  force,  pitching  moment  and  rolling 
moment . 

A  linear  variable  differential  transducer  (LVDT)  was  mounted  between  a  crank 
on  the  shaft  of  the  excitation  mechanism  and  a  fixture  on  the  sidewall 
turntable  to  measure  both  the  mean  and  unsteady  component  of  the  relative 
position  of  the  balance  beam  with  respect  to  the  sidewall  turntable.  Mean 
incidence  of  the  turntable  was  measured  by  a  very  sensitive  accelerometer 
attached  to  the  turntable,  measuring  a  component  of  the  gravitational 
acceleration. 

The  basic  wing  panel  was  provided  with  four  chordwise  rows  (Sections  1  to  4) 
and  two  spanwise  rows  (Sections  6  and  7)  of  miniature  pressure  transducers  as 
shown  in  Figure  2.  A  third  spanwise  row  (Section  5)  was  provided  in  the  strake 
(Figure  3).  The  pressure  points  of  the  chordwise  rows  were  located  toward  the 
wing  tip  in  order  to  concentrate  instrumentation  in  the  region  where  phenomena 
like  shock-induced  separation  and  leading  edge  separation,  take  place.  The 
three  spanwise  rows  were  intended  to  observe  shock-vortex  interactions  at 
higher  incidences.  Sections  5  and  6  (40.42  and  65.88  per  cent  root  chord) 
correspond  to  spanwise  sections  in  the  model  of  Reference  1.  All  pressure 
transducers  were  mounted  such  that  they  were  electrically  isolated,  free  of 
local  model  deformations,  and  not  influenced  by  model  accelerations.  The 
sensitivity  of  the  pressure  transducers  showed  a  small  variation  with 
temperature.  By  measuring  the  model  temperature  with  two  temperature  sensors. 
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the  correct  sensitivity  of  the  pressure  transducers  could  be  selected  for 
processing  recorded  electrical  signals  into  pressure  data  (see  also  Section 
4.2.2)  . 

Accelerometers,  measuring  accelerations  in  z-direction,  were  located  at  12 
positions  on  the  basic  wing  as  shown  in  Figure  2  and  at  3  positions  on  the 
strake  as  shown  in  Figure  3 . 

The  location  of  the  instriunentation  was  included  in  the  printout  of  each  data 
point.  As  a  reference  this  information  is  summarized  in  Figure  8  and  Table  1. 
In  total,  95  pressure  transducers,  1  displacement  transducer,  2  temperature 
sensors,  15  accelerometers  and  3  balance  components  were  measured  during  the 
sis  test. 

During  the  test  some  sensors  were  not  working  properly.  Table  2  gives  an 
overview. 

2 . 3  Measurement  system 

The  wind  tunnel  test  was  performed  using  the  fully  digital  computer  controlled 
measurement  system  PHAFAO  (Processor  for  Harmonic  And  Random  Oscillations) . 

The  system  acquired  128  signals  simultaneously. 

The  measurement  system  as  used  during  the  test  is  outlined  in  Figure  9.  Three 
main  parts  can  be  distinguished: 

1  Multichannel  Conditioning  Units  (MCCU's). 

2  Multichannel  Front  Ends  (MCFE's). 

3  Computers  to  control  the  MCCU's  and  MCFE's. 

In  the  MCCU's  the  filters  for  splitting  the  signal  into  an  ac  and  a  dc  part, 
as  well  as  the  preamplifier  and  the  ac  and  dc  postamplifiers  are  integrated  in 
one  module  for  each  two  channels,  together  with  a  shared  excitation  source. 
Eight  of  these  modules  are  integrated  in  one  Multichannel  Conditioning  Unit 
(MCCU) ,  which  is  suitable  for  16  input  signals.  Each  MCCU  has  its  own  ADC  and 
computer  interface.  The  dc  measurement  is  done  in  the  MCCU's.  As  a  special 
feature  the  dc  component  of  all  channels  can  be  stored  in  a  track/hold  circuit 
and  the  high-pass  filter  can  be  bypassed.  In  this  way,  a  fixed  voltage  can  be 
subtracted  from  the  preamplified  ac+dc  signal.  The  aim  of  this  mode  is  the 
possibility  to  measure  and  amplify  slowly  varying  input  signals  with  a  large 
dc  component.  This  mode  is  called  the  subtract  mode  of  the  MCCU's;  it  was 
used  in  the  SiS  test  during  the  simulation  of  maneuvers.  The  MCCU's  were 
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Fig.  8  Location  of  instrumentation  in  outboard  wing  panel  (see  also  table  1) 


Table  la  Coordinates  of  accelerometers  in  wing  and  strake 
(measures  in  mm;  see  also  Figure  8) 


NR. 

X 

11 

-425.6 

12 

-215.6 

13 

167.4 

21 

-138.6 

22 

-  46.6 

23 

121.4 

31 

-  74.6 

32 

-  10.6 

33 

141.4 

41 

29.4 

42 

89.4 

43 

152.4 

51 

85.0 

52 

121.4 

53 

157.4 

-12.0 
-12.0 
-  12.0 

-116.9 

-116.9 

-116.9 

-189.9 

-189.9 

-189.9 

-304.9 

-304.9 

-304.9 

-374.9 

-374.9 

-374.9 


Table  lb  Coordinates  of  temperature  sensors 
(measures  in  mm;  see  also  Figure  8) 


NR. 

X 

y 

1 

96.4 

-194.9 

2 

147.4 

-371.9 
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Table  Ic  Coordinates  of  pressure  points  of  Sections  1  and  2 
(measures  in  ram;  see  also  Figure  8) 


PRESSURE  SECTION  1 

chord  =  300.65 

PRESSURE  SECTION  2 

chord  =  246.21 

y  =  -209-06 

y/b  =  0.500 

y  =  -273.97 

y/b  =  0.656 

nr 

X 

X 

up 

low 

chord 

101 

2.00 

-75.31 

-21.93 

102 

5.00 

-66.29 

-14.55 

103 

151 

10.00 

-51.25 

-2.24 

104 

15.00 

-36.22 

10.07 

105 

18.00 

-27.20 

17.46 

152 

20.00 

-21.19 

22.38 

106 

30.00 

8.88 

47.00 

107 

153 

40.00 

38.95 

71.62 

108 

50.00 

69.01 

96.24 

109 

154 

60.00 

99.08 

120.86 

110 

70.00 

129.15 

145.47 

111 

79.00 

156.21 

167.63 

155 

80.00 

159.22 

170.09 

112 

82.50 

166.73 

176.25 

113 

85.00 

174.25 

182.40 

114 

90.00 

189.28 

194.71 

115 

95.00 

204.32 

207.02 
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Table  Id  Coordinates  of  pressure  points  of  Sections  3  and  4 
(measures  in  ram;  see  also  Figure  8) 


PRESSURE  SECTION  3 

chord  =  194.13 

PRESSURE  SECTION  4 

chord  =  144.42 

y  =  -336.06 

y/b  =  0.804 

y  =  -395.32 

y/b  =  0.946 

nr. 

X 

X 

up  low 

chord 

301 

2.0 

29.13 

77.86 

302 

5.0 

34.95 

82.19 

303 

10.0 

44.65 

89.40 

304 

15.0 

54.35 

96.92 

305 

18.0 

60.18 

100.95 

306 

30.0 

83.46 

118.27 

307 

40.0 

102.87 

132.70 

308 

50.0 

122.28 

147.13 

309 

60.0 

141.68 

161.56 

310 

70.0 

161.09 

175.99 

311 

79.0 

178.55 

188.98 

312 

90.0 

199.90 

204.86 
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Table  le  Coordinates  of  pressure  points  of  Section  5 
(measures  in  mm;  see  also  Figure  8) 


PRESSURE  SECTION  5 

b  =  82.7 

X  =  -269.6 

(40.42  %  cr) 

nr. 

y 

up 

Span 

501 

6.62 

5.47 

502 

20.43 

16.90 

503 

34.05 

28.16 

504 

47.67 

39.42 

505 

54.49 

45.06 

506 

61.29 

50.69 

507 

68.10 

56.32 

508 

74.91 

61.95 

509 

81.72 

67.59 

510 

88.53 

73.21 

I 

I 
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Table  If  Coordinates  of  pressure  points  of  Section  6 
(measures  in  mm;  see  also  Figure  8) 


PRESSURE  SECTION  6 

b  =  233.73 

X  =  -  60.62 

(65.88  %  cr) 

nr. 

% 

y 

up 

span 

601 

38.90 

90.91 

602 

42.93 

100.35 

603 

46.93 

109.70 

604 

50.99 

119.17 

605 

59.03 

137.96 

606 

67.07 

156.77 

607 

71.11 

166.21 

608 

75.56 

176.21 

609 

80.00 

186.99 

610 

84.44 

197.37 

102 

89.45 

209.06 
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Table  Ig  Coordinates  of  pressure  points  of  Section  7 
(measures  in  mm;  see  also  Figure  8) 


PRESSURE  SECTION  7 

b  =  417.90 

X  =  100.71 

(85.54  %  cr) 

nr. 

y 

up 

span 

701 

22.71 

94.90 

702 

28.21 

117.90 

703 

33.72 

140.90 

704 

39.26 

164.07 

705 

44.69 

186.74 

109 

50.03 

209.06 

706 

55.28 

231.03 

707 

60.46 

252.67 

208 

65.56 

273.97 

708 

70.59 

294.99 

709 

75.54 

315.69 

307 

80.42 

336.06 

710 

85.22 

356.12 

711 

90.19 

376.90 

405 

94.60 

395.32 

Table  2  Overview  of  sensors  which  did  not  work  properly 


PRESSURE  TRANSDUCERS 

DATA  POINT 

ACCELEROMETERS 

NUMBER 

St. 

unst . 

(dpn) 

105 

105 

all  dpn's 

210 

all  dpn's 

22 

all  dpn's 

32 

all  dpn's 

21 


Rg.  9  Measurement  system  during  test  on  SiS  configuration 
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designed  by  NLR.  For  other  features  of  the  MCCU's  the  reader  is  referred  to 
Reference  5 . 

Each  Multichannel  Front-End  (MCFE)  has  a  750  kHz  ADC  and  a  crystal  controlled 
clock.  In  every  MCFE  64  modules  are  integrated.  Each  module  has  a 
computer-programmable  (antialiasing)  low-pass  filter  and  a  pre-  and  postfilter 
amplifier.  The  front-end  modules  are  also  provided  with  a  sample/hold  circuit 
to  enable  real  simultaneous  sampling  of  the  different  channels.  Special 
synchronization  between  the  two  MCFE's  was  used  to  ensure  simultaneous 
sampling  of  all  signals.  The  highest  possible  sample  rate  of  the  PHARAO 
configuration  depicted  in  Figure  9  was  about  11  kHz.  Especially  at  high 
frequencies,  the  data  flow  between  MCFE's  and  computer  would  be  too  large. 
Therefore  the  MCFE's  are  equipped  with  1  Mbyte  local  memory,  which  acts  like 
a  buffer  between  the  measured  samples  and  the  memory  of  the  host  computer.  In 
one  of  the  MCFE's  a  quadruple  DAC  was  embedded.  This  DAC  was  fully 
programmable,  providing  both  the  harmonic  signals  for  model  oscillation  and 
the  signals  which  represent  maneuver  inputs. 

The  MCFE ' s  were  operated  remotely  controlled  according  to  the  scheme  in  Figure 
9.  The  measuring  process  was  governed  by  the  Computer  of  PHARAO  (COP) .  It 
controlled  directly  the  MCFE's.  Control  of  MCCU's  was  exercised  by  the  Data 
Acquisition  computer  of  the  wind  Tunnel  (DAT) ,  but  only  on  instruction  of  the 
PHARAO  computer.  In  this  way  ranging,  setting  the  MCCU  mode,  overload 
detection,  and  ADC  scans  for  all  MCCU's  simultaneously,  could  be  performed 
independently  by  the  tunnel  computer.  On  request  by  COP,  the  data  acquired  by 
DAT  were  sent  to  COP.  At  each  data  point  the  measured  time  traces  and 
corresponding  spectra  as  calculated  by  a  microprogrammed  FFT  coprocessor  card 
in  the  COP  system,  were  stored  on  an  optical  disk.  To  this  data  were  added  the 
dc  data  and  (amplifier-)  settings  obtained  from  DAT  and  the  corrected  tunnel 
data  obtained  from  the  data  processing  computer  of  the  wind  tunnel  (DPT) . 
Processing  of  the  data  to  (quick-look)  results  was  done  by  a  separate  computer 
which  had  direct  access  to  the  data  on  the  optical  disk  of  COP. 
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3  PREPARATORY  TESTS 

Before  the  tunnel  entry  of  the  model,  the  natural  frequencies,  corresponding 
vibration  modes  and  generalized  masses  were  determined  by  the  impact  method  in 
a  separate  test.  The  model  was  mounted  on  its  main  balance  beam.  The  balance 
beam  was  mounted  on  the  shaft  which  was  supported  by  bearings  mounted  at  the 
turntable  of  the  tunnel  sidewall.  The  (very  stiff)  turntable  was  supported  in 
a  special  frame  which  provided  the  possibility  of  doing  these  tests  outside 
the  wind  tunnel . 

With  the  results  of  these  tests  it  was  possible  to  choose  oscillation 
frequencies  for  the  wind  tunnel  test  program  in  such  a  way  that  (higher 
harmonics  of)  the  driving  frequency  were  not  very  close  to  natural  frequencies 
of  the  model.  The  results  of  the  vibration  analysis  were  also  used  as  input 
for  the  flutter  calculations  which  were  done  before  the  wind  tunnel  test. 

As  an  example.  Figure  10  shows  the  frequency  response  function  of  an  impact 
and  the  response  for  the  SiS  configuration.  The  different  peaks  associated 
with  the  natural  vibration  modes  can  be  distinguished  clearly.  These  modes  are 
shown  in  Figures  lib  through  lid,  while  Figure  11a  presents  the  undeformed 
grid. 

Measuring  with  balances  in  unsteady  wind  tunnel  tests  requires  the  subtraction 
of  the  inertial  loads  from  the  total  loads  to  obtain  the  aerodynamic  loads. 

For  harmonic  excitation  with  frequencies  below  8  Hz,  the  model  was  considered 
rigid  and  the  (small)  inertial  loads  were  easily  determined  from  known  mass 
properties  of  the  model  and  angular  accelerations  determined  from  the 
displacement  transducer  output.  For  the  processing  of  data  of  runs  with 
harmonic  excitation  at  frequencies  above  8  Hz  or  maneuver  data  the  "mode  fit" 
procedure  (see  Section  4. 2. 3.1)  was  used.  For  this  procedure  it  was  necessary 
to  determine  the  relation  between  the  output  of  all  accelerometers  and  the 
main  balance  for  all  natural  frequencies.  These  relations  were  measured 
during  a  preparatory  test  in  which  the  model  was  oscillated  harmonically  at 
all  natural  frequencies  below  250  Hz.  Also  these  relations  were  measured  (at 
10  Hz)  for  a  pure  pitch  motion  of  the  rigid  model. 

Note  that  the  deformations  of  the  balance  are  presented  for  each  data  point 
(see  e.g.  Table  6).  For  this  purpose  the  stiffness  matrix  of  the  balance  was 
determined  in  a  separate  preparatory  jig  test,  enabling  the  calculation  of  the 
deformations  during  the  data  reduction. 
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Fig.  10  Frequency  response  function  of  an  impact  at  point  84  and  the  response 
at  point  1 13  (see  Fig.  1  la) 
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FRONT  VIEW 


Fig.  1 1a  Undeformed  grid  used  for  impact  testing 
•  point  84,  X  point  113 
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Fig.  1 1b  Natural  vibration  mode  at  91.2  Hz:  torsion  of  the  main  balance  and  pitch 
rotation  of  the  wing  with  (small)  additional  bending  of  wing  tip  and  strake; 
(see  indicated  nodal  line  (---)) 
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Rg.  11c  Natural  vibration  mode  at  136.3  Hz:  bending  combined  with  torsion,  of  the 
wing  tip  and  a  pitching  motion  of  the  root  section  with  opposite  phase  as 
the  tip  rotation  (see  indicated  nodal  line  (---)) 
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Fig.  1 1d  Natural  vibration  mode  at  166.5  Hz:  second  bet 
outboard  wing  panel 
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of  the 


4  PROCEDURES 


4 ■ 1  Data  acquisition 

During  harmonic  oscillation,  signals  were  sampled  with  a  frequency  being  a 
multiple  of  the  model  frequency.  Both  frequencies  were  perfectly  synchronized 
as  they  were  generated  by  the  same  clock  (see  Section  2.3). 

The  cut-off  frequency  of  the  first  group  of  64  antialiasing,  equal  time  delay 
(ETD) ,  60  dB/octave  filters  which  were  programmable  in  steps  of  1  Hz,  was 
always  50  %  of  the  bandwidth.  The  second  group  of  64  antialiasing,  elliptic 
(ELL),  80dB/octave  filters  had  steeper  filter  characteristics. 

Therefore,  the  largest  cut-off  frequency  below  80  %  bandwidth  was  chosen.  Due 
to  the  large  increments  (100  Hz)  in  the  programmable  cut-off  frequency,  this 
frequency  was  not  always  close  to  80%  bandwidth.  Table  3  presents  an  overview 
of  the  different  combinations  of  the  model  frequency,  sample  frequency, 
bandwidth  and  filter  setting  used  for  the  data  acquisition. 

For  all  data  points,  including  the  maneuver  simulations,  always  4096  samples 
were  stored  in  the  measurement  file  for  each  measurement  channel.  Before  and 
after  each  series  of  measurements  the  output  of  the  sensors  at  zero-wind 
condition  was  stored  in  a  so-called  zero-measurement  file. 

4 ■ 2  Data  processing 

4.2.1  General 

For  the  data  points  with  harmonic  oscillation  of  the  model,  the  prime  interest 
was  in  the  different  harmonic  components  (multiples  of  the  model  frequency) . 
Therefore,  Phase  Locked  Time  Domain  Averaging  (PLTDA)  was  applied  on  the  time 
traces  stored  in  the  measurement  files.  Further  processing  was  done  with  the 
averaged  time  trace  of  one  cycle  of  the  model  motion.  After  the  correction  for 
the  characteristics  of  the  antialiasing  filters  and  the  determination  of  the 
squared  terms  in  the  time  domain,  Fourier  transformation  to  the  frequency 
domain  was  applied.  Application  of  sensor  characteristics,  corrections  for 
zero  measurements  and  determination  of  inertia  loads  on  the  balance  were 
performed  in  the  frequency  domain.  For  data  points  with  the  model  oscillating 
harmonically,  the  result  was  the  mean  value,  the  first  harmonic  and  the  second 
to  the  sixty-fourth  harmonic  component  in  the  spectrum  (realize  that  128 
samples  per  cycle  were  measured) . 

For  the  processing  of  the  maneuvers,  the  procedure  was  in  principle  the  same 
as  above,  except  the  fact  that  no  averaging  was  done.  So,  one  block  of  4096 
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Table  3  Data  acquisition  parameters.  For  5.7  to  15.2  Hz  only  harmonic 

oscillations  were  executed  (see  Tab.  4  (test  program));  for  3.8  Hz 
only  maneuvers  were  done  (see  Tab.  5  (test  program)  ) ,  for  this  case 
the  "frequency"  refers  to  1/timescale  of  the  full  (1-cos)  maneuver. 


f 

m 

f 

s 

128*f 

m 

Bandwidth 

(BW) 

0.5*f 

s 

64  *f 

m 

f 

c 

ETD 

50  X  BW 

0.25*f 

s 

32  *f 

m 

f 

c 

ELL 

3.8 

486.4 

243.2 

121 

100  (  41X  BW) 

5.7 

729.6 

364.8 

182 

200  (  55X  BW) 

7.6 

972.8 

486.4 

243 

300  (  62X  BW) 

11.4 

1459.2 

729.6 

364 

500  (  69X  BW) 

15.2 

1945.6 

972.8 

486 

700  (  72X  BW) 
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samples  for  each  sensor  was  processed  as  a  whole  (realize  that  for  the 
processing  of  the  data  points  with  oscillating  model  only  128  "averaged 
samples"  representing  one  cycle  of  the  model  motion,  were  processed  for  each 
sensor) . 

For  the  data  points  with  harmonic  excitation  only  the  first  16  out  of  64 
spectral  lines  of  the  results  were  considered  to  be  relevant.  Therefore  only 
these  values  were  stored  in  files  to  be  printed  and  plotted  for  data 
presentation.  For  these  spectral  lines  the  correction  due  to  the  division  by 
the  antialiasing  filter  characteristic  was  small.  For  the  maneuvers,  in 
principle  all  spectral  lines  had  to  be  processed  to  get  the  right  time  trace 
when,  after  corrections,  the  spectrum  was  Fourier  transformed  again.  Because 
the  magnitude  of  the  filter  characteristic  was  very  small  for  the  right  hand 
side  of  the  spectrum,  resolution  of  the  measurement  became  poor  and  dividing 
by  the  small  number  (of  the  filter  characteristic)  would  yield  a  noisy  result. 
Therefore  the  spectral  lines  above  1.5  times  the  cut-off  frequency  of  the 
elliptic  filters  (1.5*100=150  Hz)  were  not  used  to  reconstruct  the  corrected 
time  traces.  Also  for  the  data,  acquired  with  the  use  of  the  equal  time  delay 
filters,  spectral  lines  above  150  Hz  were  neglected  during  the  maneuver 
processing. 

4.2.2  Pressures 

Before  the  test  the  transducers  were  calibrated  at  three  different 
temperatures  within  the  range  expected  during  the  wind  tunnel  test.  Although 
the  transducers  were  so-called  temperature  compensated,  the  response  on 
temperature  steps  was  established  in  an  8-hours  test,  to  enable  correction  for 
zero  drift  due  to  temperature  effects.  Zero  drift  in  time  was  assumed  to  be 
linear  between  the  two  measurements  at  wind-off  condition  (zero-measurements), 
which  ,-.'ere  performed  before  and  after  each  series  of  measurements.  At  each 
data  point  the  temperature  of  the  pressure  transducers  was  measured  and  the 
corresponding  calibration  curves  were  used  for  processing.  During  the  tests, 
some  pressure  transducers  were  not  working  properly  (Table  2  gives  an 
overview) ;  their  presence  is  ignored  in  the  data  processing.  For  details  of 
the  pressure  integration  to  section  coefficients,  see  Section  4. 2. 3. 2. 
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4.2.3  Overall  loads 


4 . 2 . 3 . 1  Main  balance 

To  obtain  the  (unsteady)  aerodynamic  loads,  the  inertia  loads  must  be 
subtracted  from  the  total  loads  measured  by  the  balance.  As  the  natural 
frequencies  of  the  wind  tunnel  model  were  within  the  range  of  planned 
frequencies  (including  higher  harmonics  of  those  test  frequencies)  to  be 
investigated,  the  model  could  not  be  considered  rigid  anymore.  Therefore, 
inertia  loads  were  determined  with  a  procedure  that  accounts  for  the 
flexibility  of  the  wing  and  main  balance.  This  procedure  was  called  the  "mode 
fit  procedure" .  It  is  based  on  the  assumption  that  any  motion  of  the 
harmonically  excited  wing  model  can  be  expressed  as  the  weighted  sum  of  the 
driven  pitching  mode  and  a  limited  number  of  natural  vibration  modes. 

During  preparatory  tests  the  relation  between  balance  loads  and  natural 
vibration  modes  as  measured  by  the  accelerometers  in  the  model  was  established 
by  exciting  the  model  at  all  its  natural  modes  with  frequencies  below  250  Hz 
(see  Section  3) .  Also  the  balance  loads  at  a  low  pitching  frequency  (10  Hz) 
were  determined. 

During  the  wind  tunnel  test  the  actual  vibration  mode  was  fitted  as  a  linear 
combination  of  the  natural  vibration  modes.  As  the  balance  loads  for  each 
natural  mode  were  known  (see  Section  3),  the  total  inertia  loads  correction 
was  determined  as  a  weighted  linear  summation  of  those  loads. 

Preparatory  tests  were  performed  to  check  the  applicability  of  this  method  for 
the  type  of  "beam  balance"  as  used  in  this  test.  These  tests  were  reported  in 
References  6  and  7 . 

For  harmonic  excitation  with  frequencies  below  8  Hz  the  model  was  considered 
as  rigid  and  (small)  inertia  loads  were  easily  determined  from  known  mass 
properties  of  the  model  and  angular  accelerations  determined  from  the 
displacement  transducer  output.  Higher  harmonics  of  those  data  points  were 
processed  in  the  same  way.  Above  8  Hz  the  excitation  frequency  or  higher 
harmonics  of  the  excitation  frequency  approached  the  natural  frequencies  of 
the  model.  Therefore  all  harmonics  at  tests  with  excitation  frequencies  higher 
than  8  Hz  were  processed  with  the  "mode  fit  procedure". 

Although  the  time  scale  of  the  maneuvers  was  fairly  large  (the  length  of  a 
(1-cos)  input  corresponded  to  a  frequency  of  3.8  Hz),  all  spectral  lines  up  to 
150  Hz,  which  was  well  beyond  the  first  natural  frequency  of  the  model,  were 
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processed.  To  be  consistent,  for  all  spectral  lines  the  inertia  loads  were 
determined  with  the  "mode  fit  procedure" . 

4  ■  2  ■  3  ■  2  Section  coefficients 

By  integrating  (trapezoidal  rule)  the  measured  steady  and  unsteady  pressures, 
section  coefficients  were  obtained.  For  the  chordwise  sections  between  the 
leading  edge  and  the  first  pressure  transducer  a  constant  static  and  unsteady 
pressure  was  assumed.  At  the  trailing  edge  the  static  pressure  coefficient  was 
assumed  to  be  zero.  Between  the  trailing  edge  and  the  last  pressure  transducer 
a  constant  unsteady  pressure  was  assumed.  In  the  same  way  spanwise  sections 
were  integrated  from  the  plane  of  symmetry  to  the  tip.  Between  the  symmetry 
plane  and  the  first  pressure  transducer  a  constant  steady  and  unsteady 
pressure  was  assumed.  At  the  tip  the  static  pressure  coefficient  was  assumed 
to  be  zero.  Between  the  tip  and  the  last  pressure  transducer  a  constant 
unsteady  pressure  was  assumed.  The  result  of  the  pressure  integration  was  NOT 
multiplied  by  1/tu  or  2/:c. 

4.2.4  Tunnel  corrections 

For  the  test  on  the  simple  straked  delta  wing  the  maximum  value  of  the 
blockage  and  upwash  were  estimated  to  be  3  %  of  the  dynamic  pressure  and  0.5 
deg.  As  tests  on  models  like  this  simple  straked  delta  wing  at  such  extreme 
test  conditions  are  rare,  insufficient  reference  to  define  blockage  and  upwash 
as  a  function  of  other  parameters  like  frontal  area  and  Cd  exists.  Therefore, 
no  blockage  or  upwash  corrections  were  applied. 


5  TEST  PROGRAM 


The  sis  model  was  tested  at  only  three  Mach  numbers: 

M  =  0.225  To  establish  continuity  with  the  low  speed  test  of  Reference  1. 

M  =  0.600  As  an  intermediate  Mach  number  to  study  the  change  in  vortex 

flows  with  increasing  compressibility  effects. 

M  =  0.900  To  study  the  shock/vortex  interactions. 

All  three  Mach  numbers  were  tested  at  a  Reynolds  number  (based  on  the  root 
chord)  of  8.0*10®.  Table  4  presents  the  program  of  the  test  with  the 
oscillating  model.  At  M=  0.225  also  the  tunnel  was  evacuated  to  obtain  the 
Reynolds  n\amber  (based  on  the  root  chord)  of  Reference  1,  namely  3.8*10®.  At 
M=  0.900  a  series  of  incidences  was  also  tested  at  a  Reynolds  number  of 
14*10®.  The  same  (reduced-) frequency-amplitude  combinations  were  used  as  in 
the  test  described  in  Reference  1. 

The  tests  on  the  SiS  configuration  were  started  with  two  incidence  sweeps  at 
the  Mach/Reynolds  conditions  of  Reference  1,  with  and  without  the  filler  plate 
(see  Section  2.2.3).  The  results  with  the  filler  plate  installed  showed  the 
best  agreement  with  the  measurements  on  the  full-span  model  with  the  same 
planform  reported  in  Reference  1.  It  was  decided  therefore  to  perform  all  the 
other  measurements  (dpn's  >  60)  with  the  filler  plate  installed. 

Also  maneuver  sequences  were  performed  with  time  traces  of  the  model  incidence 
as  presented  in  Figure  12 .  Table  5  presents  an  overview  of  the  conditions  for 
which  maneuver  simulations  were  done;  in  this  table  alpha  does  NOT  indicate 
the  starting  angle  of  attack  but  the  mean  value  of  the  adjusted  minimum  and 
maximum  angle  of  attack  during  the  sequence,  while  dalpha  indicates  half  the 
difference  between  minimum  and  maximum. 

In  the  test  programs  adjusted  values  have  been  collected.  Measured  values  of 
a,  da,  etc.  are  included  in  the  data  presentation. 
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Table  4a  Test  program,  Mach  =  0.225,  harmonic  oscillations 


#  without  filler  plate 


Table  4b  Test  program,  Mach  =  0.600,  harmonic  oscillations 


Table  4c  Test  program,  Mach  =  0.900,  harmonic  oscillations 


Mach 

0.900 

Re*10-‘ 

-  8.0 

-  14.0 

£req. 

5.7 

7.6 

11.4 

15.2 

5.7 

dalpha 

o 

i/i 

4.0 

8.0 

4.0 

8.0 

2.0 

2.0 

in 

o 

alpha 

4.0 

499 

527 

5.0 

500 

528 

6.0 

501 

566 

574 

584 

592 

601 

609 

617 

529 

7.0 

502 

530 

8.0 

503 

531 

9.0 

504 

533 

10.0 

505 

567 

575 

585 

593 

602 

610 

618 

534 

10.5 

535 

11.0 

506 

536 

11.5 

537 

12.0 

507 

538 

12.5 

539 

13.0 

508 

540 

14.0 

509 

568 

576 

586 

594  595 

603 

611 

619 

541 

15.0 

510 

542 

16.0 

511 

543 

17.0 

512 

544 

18.0 

513 

569 

579 

587 

596 

604 

612 

620 

545 

19.0 

514 

548 

20.0 

515 

547 

21.0 

516 

549 

22.0 

517 

570 

580 

588 

597 

605 

613 

621 

526  550 

23.0 

518 

551 

24.0 

519 

565 

583 

591 

600 

608 

616 

25.0 

520 

553 

26.0 

521 

571 

581 

589 

598 

606 

614 

622 

554 

27.0 

522 

555 

28.0 

523 

556 

29.0 

524 

557 

30.0 

525 

572 

582 

590 

599 

607 

615 

623 

558 

32.0 

559 

34.0 

560 

36.0 

561 

38.0 

562 

40.0 

563 

42.0 

44.0 

46.0 

48.0 

564 
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Table  5  Test  program,  maneuvers;  "freq"  =  3.8  Hz 


Mach 

0.225 

0. 

600 

0.900 

Re*10‘‘ 

-  3.8 

• 

8.0 

8.0 

-  8.0 

dalpha 

8.0 

8.0 

16.0 

8 

.0 

16.0 

8.0 

16.0 

alpha 

6.0 

235  236 

454 

455 

625  639 

237  238 

456 

457 

640  641 

239 

458 

642  643 

10.0 

240  241 

242  243 

244 

14.0 

245  246 

296  297 

459 

460 

487  488 

644  645 

662  663 

247  248 

298  299 

461 

462 

489  490 

646  647 

664  665 

249 

300 

463 

464 

491 

666  667 

18.0 

64  65 

250  251 

301  302 

67  75 

252  253 

303  304 

76  77 

254 

305 

22.0 

255  256 

306  307 

449 

450 

480  481 

624  626 

656  657 

257  258 

308  309 

451 

452 

482  483 

627  628 

658  659 

310 

453 

484 

629  630 

631  632 

633  634 

635  636 

637  638 

660  661 

26.0 

78  81 

229  230 

291  292 

82  87 

231  232 

293  294 

88 

233  234 

260  261 

262  263 

264  265 

295 

30.0 

266  267 

311  312 

465 

466 

492  493 

648  649 

268  269 

313  314 

467 

468 

494  495 

650  651 

315 

469 

470 

496 

652  653 

471 

472 

654 

473 

474 

34.0 

89  90 

276  277 

316  317 

92  93 

278  279 

318  319 

95 

280 

320 

38.0 

281  282 

475 

476 

283  284 

477 

478 

285 

479 

42.0 

286  287 

288  289 

290 
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6  PRESENTATION  OF  RESULTS 

6 ■ 1  Wind  tunnel  test  data 


Part  2  presents  five  harmonics  of  five  selected  (see  Section  7)  data  points  of 
the  test  in  plotted  and  tabulated  form.  Table  6  and  Figure  13  show  the  format 
in  which  the  data  are  presented. 

Each  table  presents  the  test  conditions  (alpha, Mach, Re, Q, dalpha, frequency)  and 
one  harmonic  (>0)  component  and  the  mean  (0-th  harmonic  component)  value  of  all 
measured  quantities  for  one  data  point.  Pressures  and  section  loads  are 
presented  in  nondimensionalized  form.  The  inertia  loads  were  subtracted  from  the 
balance  loads,  before  overall  aerodynamic  loads  were  non-  dimensionalized  to 
aerodynamic  coefficients.  Also  the  ratio  between  the  amplitude  of  the  unsteady 
aerodynamic  loads  and  the  amplitude  of  the  inertia  loads  is  indicated  (when  **** 
are  printed  the  ratio  is  larger  than  10,000).  The  angular  deflections, 
calculated  from  the  total  balance  loads  and  the  stiffness  matrices,  are 
presented  to  the  right  of  the  balance  loads.  To  the  right  of  the  accelerations 
the  vibration  mode  of  each  section  equipped  with  accelerometers  is  expressed  in 
heave  at  the  pitching  axis  and  rotation  about  this  axis.  Each  plot  presents  one 
harmonic  (>0)  component  of  the  pressure  distribution;  the  mean  (0-th  harmonic 
component)  pressure  distribution  is  always  included  as  a  reference.  The  data  of 
the  five  harmonics  of  the  five  selected  data  points  are  also  available  on  a 
floppy  disk. 

Part  2  also  presents  plots  of  overall  quantities  like  normal  force  and  pitching 
moment  coefficients  versus  incidence,  in  which  the  results  of  a  large  number  of 
data  points  have  been  collected,  including  the  five  selected  data  points. 

Part  3  presents  the  response  on  a  maneuver  sequence  in  plotted  form  for  three 
selected  (see  Section  7)  data  points.  All  the  quantities  which  were  presented 
for  the  data  points  with  harmonically  oscillating  model,  together  with  inertia 
loads  as  separate  quantities,  are  presented  both  as  a  time  trace  and  as  a 
cross-plot  with  the  incidence  signal.  As  an  example.  Figure  14  shows  the 
response  of  the  main  balance  in  terms  of  overall  load  coefficients  and  angular 
deflections  for  the  four  parts  of  a  maneuver  sequence.  The  data  presented  in  the 
plots  are  also  made  available  as  ASCII  files  on  floppy  disks. 

When  a  sensor  did  not  operate  properly,  no  value  is  presented  in  the  plots, 
prints  and  ASCII  files.  Readers  who  are  interested  to  obtain  part  or  all  of  the 
measured  data  should  contact  Lockheed  Forth  Worth  Company  or  NLR. 
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Table  6  Example  of  (quick  look-)  printout  of  the  mean  and  first  harmonic 
component  of  the  data  measured  at  one  data  point  (page  1/3) 


Unsteady  Transonic  Delta  Program 


test  conditions 

Simple  Strahe  coniiguxation 

alpha  -  21.914  dag 

Q  -3.138  kPa 

Mach  -  .223 

Ptot  =  92.051  kPa 

Re*10-‘  =.  3.804 

Itot  =•  288.834  K 

dalpha  -  .546  dag 

fraq  -  5.700  Hz 

k  -  . 193 

harm  -  1 

BALANCE  LOADS 

aerodynamic  or  <total>  coafficiants 

aaro 

angular  deflections  Cdag] 

position 

cotap. 

Zero 

Ra  1 

Iffl  1 

inertia 

til 

Zaro 

Ra  1 

Im  1 

main 

at 

Cm 

Cl 

1.24332 

.07572 

-.43299 

2.30012 

.32252 

-.52972 

.  70846 
-.03933 
-.35483 

1197.71 

103.80 

327.90 

-.026 

-.034 

-.002 

.000 

.000 

.000 

Accn.ti 

{AIIONS 

vibration  mode 

Phase  angle 

y 

Amplitude 

ral.  to  LTD! 

saction 

y/b 

heave  at  p.a 

pitch 

Cnml 

[m/s*2] 

Cdeg] 

ti] 

[cm) 

tdag] 

11 

-425.6 

-12.0 

5.091 

1.996 

12 

-215 . 6 

-12.0 

2.619 

2.024 

1 

2.878 

.043 

.530 

13 

167.4 

-12.0 

1.957 

-176.855 

21 

-138.6 

-116.9 

1.633 

15.502 

22 

-46.6 

-116.9 

2 

28.034 

.077 

.558 

23 

121.4 

-116.9 

1.614 

-165.610 

31 

-74.6 

-189.9 

.731 

18.711 

32 

-10.6 

-189.9 

3 

45.540 

.112 

.519 

33 

141.4 

-189.9 

1.779 

-166.231 

41 

29.4 

-304.9 

.252 

-175.904 

42 

89.4 

-304.9 

1.179 

-178.644 

4 

73.118 

.014 

.570 

43 

152.4 

-304.9 

1.815 

-178.404 

51 

85.0 

-374.9 

.840 

-172.336 

52 

121.4 

-374.9 

1.496 

-163.903 

5 

89.904 

.360 

.  703 

S3 

157.4 

-374.9 

1.970 

-165.651 
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PRESSURES  section  3 

c  -  194.13  an 

y  ”-336.06  on 

ReCp  1  ImCp  1 


PRESSURES  section  ^ 

c  -  144.42  1 

y  --395.32  1 

RaCp  1  ImCp  1 


404 

15.00 

*05 

18.00 

406 

30.00 

*07 

40.00 

408 

50.00 

*09 

60.00 

*10 

70.00 

*11 

79.00 

*12 

90.00 

1.49* 

.798 

-.020 

-.376 

-.707 

-1.079 

-.978 

-.491 
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Table  6  Example  of  (quick  look-)  printout  of  the  mean  and  first  harmonic 
component  of  the  data  measured  at  one  data  point  (page  3/3) 

Unsteady  Transonic  Dslta  Program  npH  >  g 


FKtSSUEtS  section  5 


b  “  82.70  ins 

X  "“269.60  on 


ReCp  1  loCp  1 


PHJESSURES  section  6 


-1.680 

-4.728 

-1.906 

-5.123 

-1.900 

-5.797 

-1.73*  • 

-6.660 

-1.329 

-7.003 

-1.407 

-8.606 

-1.809 

-9.6*6 

-2.844 

-5.587 

-3.437 

4.132 

-2.323 

2.003 

-2.012 

2.189 

233.73 

-60.62 


RaCp  1  ImCp  1 


.133 

-.173 

-.6*0 

-.382 

-.62* 

-.872 

-1.360 

-2.097 

-2.402 

-2.726 

-2.2*5 


HiESSURES  section  7 


b  -  *17.90  on 
X  "  100.71  am 


SECTION  COEITICIENTS 


-.01* 

-.308 

-1.060 

-1.518 

-l.*12 

-1.173 

-1.081 

-.938 

-.672 

-.577 

-.549 

-.495 

-.466 

-.436 

-.418 


ReCp  1 

IsCp  1 

-.061 

-.143 

.24* 

.313 

-.83* 

1.318 

-6.458 

1.368 

-11.004 

.90* 

-8.303 

-.16* 

-1.670 

-2.100 

7.98* 

-5.030 

8.206 

-4.37* 

-.220 

-.17* 

.450 

-.530 

-.294 

-.420 

-.682 

.063 

-1.86* 

1.188 

-2.110 

1.494 

Fig.  13  Example  of  (quick  look-)  plot  of  the  mean  and  first  harmonic  component  of 
the  data  measured  at  one  datapoint 
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Fig.  14  Ex3niple  of  the  response  of  the  mein  bsisnce  on  3  m3noeuver  sequence  (pege  1/4) 
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Fig.  14  Example  of  the  response  of  the  main  balance  on  a  manoeuver  sequence  (page  2/4) 
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Fig.  14  Example  of  the  response  of  the  main  balance  on  a  manoeuver  sequence  (page  3/4) 
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Fig.  14  Example  of  the  response  of  the  main  balance  on  a  manoeuver  sequence  (page  4/4) 
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6 ■ 2  Model  geometry 


Figures  1  to  7  and  Appendices  A  and  B  provide  the  most  important  geometric 
data  of  the  model.  Readers  who  are  interested  to  obtain  the  geometric  data 
(e.g.,  in  the  CATIA  or  VDAFS  format),  should  contact  Lockheed  Forth  Worth 
Company  or  NLR . 


7  BRIEF  DISCUSSION  OF  SOME  RESULTS;  SELECTION  OF  PRESENTED  DATA  POINTS 

The  rationale  for  selecting  the  data  points  presented  in  Parts  2  and  3  is 
discussed  in  this  section. 

Part  2  data  points  were  taken  from  the  oscillatory  testing  of  the  model. 

The  points  selected  for  this  part  were  chosen  to  highlight  (1)  vortex  flow 
breakdown,  and  (2)  onset  of  Shock  Induced  Trailing  Edge  Separation  (SITES)  and 
leading  edge  separation  at  transonic  speeds.  The  effect  of  Mach  number  is  the 
primary  point  of  interest  in  selecting  these  data  points. 

At  M=0.225,  data  point  151  shows  the  effect  of  the  model  oscillating  between 
14  deg  incidence  and  30  deg  incidence  at  5.7  Hz.  With  vortex  bursting  starting 
at  about  22  deg,  this  oscillation  provides  maximum  pitch  rate  at  burst  point. 
Similar  data  are  given  for  M=0.60  in  data  point  375  where  vortex  breakdown 
apparently  begins  between  23  and  24  deg.  The  data  point  chosen  for  M=0.90  is 
605  at  the  higher  frequency  of  7.6  Hz  which  provides  an  approximately  constant 
reduced  frequency  when  compared  with  data  point  375  at  Mach=0.60.  In  the  case 
of  data  point  605,  vortex  bursting  begins  at  about  18  deg  incidence  for 
M=0.90.  As  a  reference,  the  mean  and  first  harmonic  component  of  the  pressure 
distributions  of  data  points  151,375  and  605  are  included  in  this  part  of  the 
report  (Figures  15,  16  and  17). 

The  onset  of  SITES  and  leading  edge  separation  occurs  between  10  and  12  deg  on 
the  simple  straked  wing  at  Mach=0.90.  In  order  to  highlight  these  transitions 
data  point  358  was  chosen  to  show  how  aerodynamics  responded  to  oscillations 
of  4  deg  amplitude  at  10  deg  mean  angle  and  M=0.60,  where  no  such  transitions 
occur.  Data  point  593  was  then  chosen  to  show  SITES/leading  edge  separation 
transition  at  M=0.90  at  the  same  conditions.  Frequency  for  the  M=0.60  case  was 
5.7  Hz  and  for  the  M=0.90  case  was  7.6  Hz  in  order  to  maintain  an 
approximately  constant  reduced  frequency.  Figures  18  and  19  show  the  mean  and 
first  harmonic  component  of  the  pressure  distributions  of  data  points  358  and 
593. 

Part  3  data  points  are  taken  from  the  large  amplitude  maneuvering  part  of  the 
simple  straked  wing  test.  The  large  amplitude  motions  of  15  deg  amplitude 
centered  on  a  mean  angle  of  22  deg  provide  a  dynamic  variation  of  flow  fields 
covering  attached,  vortex,  burst  vortex  and  developing  separated  flows  for 
incidences  from  7  to  37  deg.  The  three  Mach  numbers  are  represented  by  data 
point  306  (M=0.225),  480  (M=0.60),  and  656  (M=0.90).  In  all  cases  the 
frequency  was  held  constant  at  3.8  Hz  in  order  to  simulate  the  same  maneuver 
at  different  speeds.  Since  these  data  points  are  for  transient  and  not  for 
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Mean  and  first  harmonic  component  of  the  pressure  distribution  measured  at 
datapoint  151  (first  datapoint  selected  for  extensive  presentation  in  part  II) 


Fig.  16  Mean  and  first  harmonic  component  of  the  pressure  distribution  measured  at 
datapoint  375  (second  datapoint  selected  for  extensive  presentation  in  part  II) 
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Fig.  1 7  Mean  and  first  harmonic  component  of  the  pressure  distribution  measured  at 
datapoint  605  (third  datapoint  selected  for  extensive  presentation  in  part  II) 
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Fig.  18  Mean  and  first  harmonic  component  of  the  pressure  distribution  measured  at 
datapoint  358  (fourth  datapoint  selected  for  extensive  presentation  in  part  II) 
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Fig.  19  Mean  and  first  harmonic  component  of  the  pressure  distribution  measured  at 
datapoint  593  (fifth  datapoint  selected  for  extensive  presentation  in  part  II) 
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oscillatory  motions,  they  are  represented  in  a  time  history  format  and  thus  do 
not  have  the  harmonic  form  that  is  used  in  Part  2 .  Because  of  the  large 
quantity  of  data  within  these  data  points,  only  three  were  chosen  for  Part  3. 
In  this  part,  as  a  reference,  only  the  time  histories  of  the  balance 
quantities  are  presented  for  the  three  data  points  in  Figures  20,  21  and  22. 
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Rg.  20  Time  histories  of  the  balance  quantities  measured  at  datapoint  306 
(first  datapoint  selected  for  extensive  presentation  in  part  III)  page  1/4 
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Fig.  20  Time  histories  of  the  balance  quantities  measured  at  datapoint  306 
(first  datapoint  selected  for  extensive  presentation  in  part  III)  page  2/4 
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Ftg.  20  Time  histories  of  the  balance  quantities  measured  at  datapoint  306 
(first  datapoint  selected  for  extensive  presentation  in  part  III)  page  3/4 
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Fig.  20  Time  histories  of  the  balance  quantities  measured  at  datapoint  306 

(first  datapoint  selected  for  extensive  presentation  in  part  III)  page  4/4 
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^9-  21  Time  histories  of  the  balance  quantities  measured  at  datapoint  480 

(second  datapoint  selected  for  extensive  presentation  in  part  III)  page  1/4 
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Ftg.  21  Time  histories  of  the  balance  quantities  measured  at  datapoint  480 

(second  datapoint  selected  for  extensive  presentation  in  part  III)  page  2/4 
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Fig.  21  Time  histories  of  the  balance  quantities  measured  at  datapoint  480 

(second  datapoint  selected  for  extensive  presentation  in  part  III)  page  3/4 
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Fig.  22  Time  histories  of  the  balance  quantities  measured  at  datapoint  656 

(third  datapoint  selected  for  extensive  presentation  in  part  III)  page  1/4 
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Fig.  22  Time  histories  of  the  balance  quantities  measured  at  datapoint  656 

(third  datapoint  selected  for  extensive  presentation  in  part  III)  page  2/4 
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Fig.  22  Time  histories  of  the  balance  quantities  measured  at  datapoint  656 

(third  datapoint  selected  for  extensive  presentation  in  pari  ill)  page  3/4 
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Fig.  22  Time  histories  of  the  balance  quantities  measured  at  datapoint  656 

(third  datapoint  selected  for  extensive  presentation  in  part  III)  page  4/4 
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8  FINAL  COMMENTS 


8 ■ 1  Test  time  per  data  point 

Data  acquisition  time  was  generally  in  the  order  of  seconds.  Directly  after 
data  acquisition  the  wind  tunnel  turntable  was  rotated  to  the  next  mean 
incidence,  while  the  model  was  kept  oscillating  at  a  fixed  frequency- 
amplitude  combination.  Data  processing  was  done  when  no  acquisition  was 
needed.  For  each  data  point  quick-look  prints  and  plots  with  the  mean  and 
first  harmonic  component  of  all  sensors  were  presented  about  1  MINUTE  after 
data  acquisition  when  no  interference  with  data  acquisition  of  the  next  data 
point  occurred.  The  quick-look  format  was  the  same  as  the  format  used  for  the 
final  data  presentation. 

8 . 2  Quasi-steady  results 

Also  measurements  were  performed  at  very  low  frequencies .  No  really  steady 
measurements  were  performed;  as  a  consequence,  no  derivatives  could  be 
calculated. 

8 . 3  Accuracy 

Concernir  :he  accuracy  of  measured  aerodynamic  data,  it  is  generally 
recognized  that  it  is  extremely  difficult  to  quantify  the  accuracies  in  a 
justified  way.  Based  on  the  accuracies  of  the  instrumentation,  hardware 
components  and  the  methods  used  in  the  software,  together  with  the  results  of 
repeatability  runs  and  experience  from  previous  wind  tunnel  tests,  the 
accuracies  for  the  test  were  estimated.  The  estimates  for  the  accuracy  of  the 
mean  flow,  model  motion,  pressure  coefficients,  and  the  total  (aerodynamic  + 
inertia)  loads  on  the  balance  are  collected  in  Table  7.  The  accu  ';cy  of 
sectional  coefficients  is  determined  by  the  accuracy  of  the  pressure 
coefficients  (as  indicated  in  Table  7),  the  niomber  and  location  of  the 
pressure  transducers  and  the  gradients  in  the  (unsteady)  pressure 
distribution.  Therefore  no  single  value,  applicable  to  all  data  points,  could 
be  included  in  Table  7 .  The  accuracy  of  the  unsteady  aerodynamic  loads  as 
measured  by  the  balance  is  determined  by  the  accuracy  of  the  balance  ,  the 
ratio  between  aerodynamic  and  inertia  loads  (which  is  always  presented  in  the 
printed  tables,  see  also  Section  6.1)  and  the  accuracy  of  the  inertia  loads. 
The  accuracy  of  the  inertia  loads  on  the  balance  as  determined  by  the  "mode 
fit  procedure"  (see  Section  4. 2. 3.1)  is  determined  by  both  the  accuracy  of  the 


Table  7  Estimated  accuracies 


1.  Mean  test  conditions 


Mach  number 
Reynolds  number 
Dynamic  pressure 


+/-  0.001 
+/-  0.1 
+/-  0.2 


per  cent 
per  cent 


2.  Model  position/motion 

Turntable  position 
Steady  incidence  of  shaft 
relative  to  turntable 
Pitch  amplitude  shaft 
Pitch  deflection  of 
main  balance  beam 
Reduced  frequency 
Amplitude  of  accelerometers 


+/- 

0.002  +  0 . 0004*alpha(in  deg) 

deg 

+/- 

0.005 

deg 

+/- 

0.005 

deg 

+/- 

0.001 

deg 

+/- 

0.0005 

+/- 

(0.01*1 ampl. 1  +0.15) 

m/s 

3 .  Balance  loads 


+/•  1.0 


per  cent 


4.  Pressure  coefficients 

Steady  (mean)  +/-  0.3  per  cent 

Unsteady  +/-  0.5  per  cent 
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measured  accelerations  and  the  accuracy  of  the  procedure  which  fits  the 
measured  vibration  mode  with  a  weighted,  linear  combination  of  natural 
vibration  modes.  For  all  data  points,  processed  with  the  "mode  fit  procedure", 
the  measured  accelerations  were  fitted  with  a  mean  accuracy  better  than  1  per 
cent . 

8  ■  4  Oscillation  amplitude 

The  dalpha,  as  determined  from  the  signal  of  the  LVDT  connected  to  the  shaft 
of  the  excitation  mechanism,  is  presented  in  the  test  conditions  of  each  data 
point  (see  e.g..  Table  6,  first  page).  All  the  unsteady  signals  were 
nondimensionalized  with  this  value.  When  the  small  amplitude  of  the  torsional 
deformation  of  the  balance  beam  (see  e.g..  Table  6,  first  page)  is  added  to 
this  value,  the  amplitude  of  the  wing  pitch  motion  at  the  attachment  of  the 
wing  to  the  balance  beam  is  obtained.  When  no  torsional  deformation  of  the 
wing  is  assumed,  this  value  should  be  the  same  as  the  pitch  amplitudes  as 
obtained  from  the  accelerometers  in  the  wing  (see  e.g.,  Table  6,  lower  part 
first  page) . 

The  maximum  pitch  rate  during  the  LCO  test  (see  Section  2.2.1  and  Figure  1) 
was  much  higher  than  at  the  SiS  test.  Because  the  outboard  wing  panel  was  used 
in  both  tests,  the  range  of  the  accelerometers  (1000  m/s^)  had  to  be  high 
enough  for  the  LCO  test.  As  a  result  for  many  data  points,  especially  at  the 
SiS  test,  the  accelerometers  operated  in  a  very  small  part  of  their  range. 
Consequently,  the  pitch  amplitude  deduced  from  the  output  of  the 
accelerometers  is  not  very  accurate  (see  e.g..  Table  6,  lower  part  first  page) 
for  these  data  points. 


72 


9  CONCLUDING  REMARKS 

The  aims  mentioned  in  Section  1.2  have  been  achieved: 

A  very  large  data  base  with  steady  and  unsteady  aerodynamic  loads  on  a 
semispan  model  of  a  straked  delta  wing,  performing  both  harmonic  pitch  . 
oscillations  as  well  as  maneuvers,  was  generated.  The  data  can  be  used 
directly  for  flow  analysis  and  evaluation  of  computer  codes. 

Out  of  about  400  data  points  with  the  model  oscillating  harmonically  in 
pitch,  only  5  data  points  are  presented  extensively  in  this  report  (see 
Part  2) .  The  results  of  these  5  data  points  are  provided  in  printed  and 
plotted  form  as  well  as  in  ASCII  files  on  floppy  disk.  Out  of  about  200 
data  points  with  maneuver  simulation,  only  3  are  presented  completely  (see 
Part  3)  in  plotted  form  as  well  as  in  ASCII  files  on  floppy  disks. 
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APPENDIX  A  Axis  system  and  sign  conventions 


Figure  A.l  presents  the  body  fixed  axis  system  for  the  SiS  model: 


x-axis . 


y-axis . 


z-axis . 


In  the  wing  reference  plane  (WRP) ,  parallel  to  the  root  chordline  of 
the  basic  wing  panel  at  a  distance  of  62.3  mm.  The  WRP  is  the  plane 
defined  by  the  root  chordline  of  the  basic  wing  panel  and  the  line 
which  connects  the  0  %  chord  points  of  the  basic  wing  panel. 

In  the  wing  reference  plane,  perpendicular  with  the  x-axis,  going 
through  48.24  per  cent  of  the  root  chordline  of  the  basic  wing 
panel . 

Perpendicular  to  x-axis  and  y-axis. 


The  y-axis  coincides  with  the  rotation  axis  of  the  excitation  mechanism. 

The  z=0  plane  is  the  wing  reference  plane  (WRP) .  Both  the  root  chordline  of 
the  basic  wing  panel,  the  rotation  axis  (=  y-axis)  and  the  line  connecting  the 
0  %  chord  points  are  in  this  plane. 

The  trailing  edge  is  one  straight  line.  Due  to  the  twist,  this  line  is 
crossing  the  WRP  at  the  root  chord  of  the  basic  wing  panel.  Although  the  apex 
of  the  strake  is  in  the  WRP,  the  chordline  of  the  y=0  section  is  not  precisely 
in  the  WRP;  it  has  a  0.0803  deg  more  positive  angle  of  attack  than  the  root 
chordline  of  the  basic  wing  panel. 
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of  219.4  from  the  trailing  edge  corresponds  to  73.27%  local  chord. 


Figure  A.  1  The  body  fixed  axis  system  (dimensions  in  mm) 


APPENDIX  B  Basic  wing  panel  geometry 


This  appendix  gives  only  an  overview  of  the  geometry  of  the  wing  panel. 

The  main  dimensions  of  the  wing  panel  are  indicated  in  Figure  B.l.  The  wing 
reference  plane  (WRP)  is  defined  by  the  root  chordline  and  the  line  connecting 
the  0  %  chord  points.  The  rotation  axis  of  the  experiment  was  located  in  this 
plane.  Some  characteristics  of  the  wing  panel  are  listed  below. 

sweep  angle  :  40  deg 
of  leading  edge 

airfoil  section:  NACA  64A204  (mod.)  (all  sections);  lower  side  of  the 

sections  of  the  inboard  part  of  the  wing  is  thickened, 
twist  :  -3.0  deg;  the  y  =  -62.3  section  has  0.0  deg  incidence  with 

respect  to  the  wing  reference  plane  (WRP),  the  y  =  -417.9 
section  has  -3.0  deg  incidence;  the  panel  is  linearly  lofted 
between  root  and  tip. 
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125.36 


Figure  B.  1  Dimensions  of  the  outboard  wing  panel 

i 

I 
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